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Abstract 

New coherent states may be induced by pertinently engineering the topology 
of a network. As an example, we consider the properties of non-interacting 
bosons on a star network, which may be realized with a dilute atomic gas 
in a star-shaped deep optical lattice. The ground state is localized around 
the star center and it is macroscopically occupied below the Bose-Einstein 
condensation temperature T c . We show that T c depends only on the num- 
ber of the star arms and on the Josephson energy of the bosonic Josephson 
junctions and that the non-condensate fraction is simply given by the reduced 
temperature T/T c . 

I. INTRODUCTION 

Although it is well known that free bosons hopping on translationally invariant networks 
cannot undergo Bose-Einstein condensation at finite temperature if the space dimension d is 
less or equal to two (see Ref. [1]), very recent studies [2,3] hint to the exciting possibility that 
the network topology may act as a catalyst for inducing a finite temperature spatial Bose- 
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Einstein condensation even if d < 2. As an example of this situation we shall investigate the 
properties of non-interacting bosons hopping on a star shaped optical network, evidencing 
that - already for this very simple graph topology - one may have a macroscopic occupation 
of the ground-state at low temperatures. 

A star graph (see Fig. 1) is made of p one-dimensional chains (arms) which merge in 
one point called the center of the star. Each site i of the star arms is naturally labeled 
by two integer indices x and y where x — 0, • • • , L labels the distance from the center and 
y = 1, •••,£> labels the arms. The center of the network is denoted by O = (0, y). The 
total number of sites is N s = (pL + 1); each site on the arm is linked only to two neighbors 
whereas the center has p neighbors: thus the fact that the center has coordination number 
p is the source of spatial inhomogeneity in this lattice. In the following we shall evidence 
that bosons hopping on this graph undergo - at a certain temperature T c , which depends on 
the number of star arms - a topology induced spatial Bose-Einstein condensation in a state 
localized around the center of the star graph. 

Bosons hopping on star-shaped networks can be experimentally realized loading a di- 
lute Bose-Einstein condensate (BEC) in a suitable periodic potential, arranged to provide 
a star-like configuration. Periodic potentials are, nowadays, routinely created with two or 
more counterpropagating laser beams; one can accurately tune the height of the potential 
(which is proportional to the power of the lasers) as well as the distance between neigh- 
boring sites. When one loads an atomic BEC in a deep optical lattice, one has a bosonic 
Josephson network, i.e., an array of bosonic Josephson junctions (BJJ's). A single BJJ may 
be obtained by loading a BEC in a double well potential: the weak link between the atomic 
condensates is provided by the energy barrier between them and the dynamics of the atoms 
at T = is described by Josephson equations, obtained from the Gross-Pitaevskii equation 
[4]. Similarly, with a multi-well periodic potential, when the heights of the barriers are 
much higher than the chemical potential, the system realizes a lattice of weakly coupled 
condensates. Each bosonic Josephson junction consists of a pair of neighboring condensates: 
the tunneling rate (proportional to the Josephson energy) is easily tuned by changing the 
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power of the lasers, and it decreases if the interwell barriers increase. The properties of a 
bosonic Josephson network at T = are described by a discrete nonlinear Schrodinger equa- 
tion, obtained from the Gross-Pitaevskii equation with a periodic potential [5]. We remark 
that recently it has been showed that a 2D optical network of BEC can be described by 
the Gross-Pitaevskii equation also at finite temperature, provided that the interwell energy 
barriers and the frequency of the axial confinement are large enough [6]. 

Present-day BJJ networks are built on regular geometries: linear chains [7], squares 
[8] and cubes [9]. However, a variety of non- conventional structures may be produced by 
the standing waves of several interfering laser beams suitably placed [10]. For instance, 
a star with four arms may be realized by having two perpendicular gaussian laser beams 
superimposed to a 2D optical lattice; the center of the star should be then arranged by 
adjusting the distance between neighboring wells. 

BJJ networks may be pertinently described by the Bose-Hubbard Hamiltonian [11,12]. In 
fact, when all the relevant physical parameters are small compared to the excitation energies, 
the field operator [12], describing the condensate configuration in the BJJ network, may 
be expanded as ip(f, r) = J2j ^i( r )0i(^*) with (f>j(r) the normalized Wannier wavefunction 
localized in the j-th well and hj = a^dj the bosonic number operator. Substituting the 
expansion for Tp(r, r) in the full quantum Hamiltonian describing the bosonic system [12] 
leads to a Bose-Hubbard model 

H = -tJ2 + h.c.) + j ]Tn, - 1) . (1) 

<i,3> j 

In Eq. (1) J2<i,j> denotes a sum over all the distinct pairs of neighboring sites, aj (dj) is the 
bosonic operator which creates (destroys) a boson at site j, U = (Airh 2 'a / 'm) J drfil- (a is the 
s-wave scattering length and m is the atomic mass) and t ~ — / dr[^W<pi ■ V0j + <piV ext (f)j\, 
where 4>i an d 4>j are the Wannier functions at the neighboring sites i and j and V ext is the 
external potential confining the bosons. Since, for atomic condensates, U may be varied by 
tuning the scattering length using Feshbach resonances [1], one may assume that U <^t; for 
the sake of simplicity, in the following we shall set U — 0. 
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The filling, i.e., the average number of particles per site, is defined as / = N T /N S , where 
N T is the total number of bosons. When / 3> 1 and the fluctuations of the particle numbers 
per site are much smaller than /, one can safely substitute the operator with y/Nl e 1 ^, 
where iVj is the number of particles at the site % [14,15]. As a result, the Josephson energy 
of a single BJJ is given by 

Ej « 2tf. (2) 

It is worthwhile to remark that the Bose-Hubbard Hamiltonian (1) describes not only 
BEC's, but in general cold bosons in a deep periodic potential, provided that the temperature 
does not excite higher bands: E gap > ksT, where E gap oc ^/Vo is the band energy gap and 
V is the interwell energy barrier. In this paper we are considering a BEC in a deep optical 
lattice: thus we are also assuming that in each well one has a grain condensate, i.e., that one 
works at a temperature smaller than the temperature Tbec at which condensation occurs 
in a single well. For T < T BEC) one has a network of weakly coupled condensates, and 
the bosons, at the temperature T c defined by Eq. (14), start to macroscopically occupy the 
ground-state whose wavefunction is given by Eq. (8). We observe that for V ~ 15E R , where 
Er = h 2 /2m\ 2 with A the wavevector of the lattice beams, and for an average number of 
particles per site / ~ 200, one has Tbec ~ 500nK and T c ~ 50nK [6]. The experimental 
situation in which the total number of particles is so small that T BEC is much lower than 
T c should be also possible: in this case for T > T c one has not an array of BJJ's, but at T c 
the bosons can still macroscopically occupy the ground-state described by the eigenfunction 
(8). 

The plan of the paper is the following. In Sec. II we analyze the spectrum of bosons 
hopping on a star-shaped graph; we shall evidence that, due to the network inhomogeneity, 
an isolated localized ground-state appears in the spectrum (its effect is similar to the one 
induced by an impurity on a linear chain, see e.g. Ref. [13]); we shall also determine, as 
a function of the number of star arms p, the energy gap between the ground-state and 
the continuum states as well as the coherence length £ which characterizes the topology 
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induced spatial condensation of bosons in the center of the star graph. In Sec. Ill we 
analyze the thermodynamic properties of non-interacting bosons hopping on the star graph 
and we determine the critical temperature T c as a function of p and of the Josephson energy 
of the single BJJ; we show that the condensate fraction and the inhomogeneous spatial 
distribution of bosons over the array can be expressed in a simple way as a function of 
the scaled temperature T/T c . Section IV is devoted to some remarks on the results of our 
investigation and to a discussion of experimental settings which should enable to evidence 
the existence of the topology induced spatial BEC analyzed in this paper. 

II. SPECTRUM OF BOSONS HOPPING ON A STAR GRAPH 

The topology and geometry of a generic graph network is fully described by its adjacency 
matrix A x>mx ^ y i whose entries are 1 if (x, y; x',y') is an allowed link and otherwise; for a 
graph network, the Hamiltonian (1) is written as 

H = -t A x,y;x',y' a Xty a x f jy r. (3) 

x,y;x',y' 

The single-particle energy spectrum on a star network is found then by solving the eigenvalue 
equation [3]: 

-t A x,y; x>, y > ^e(x', y') = Eijj E (x, y), (4) 

x',y' 

with the adjacency matrix A x>y]X ^ y i given by: 

A x ,y- x ',y' — (8x',x-l + ^x',x+l) (1 — ^x,o) $y,y' + $x,0 &c',l- (5) 

We shall refer the interested reader to the Appendix for the mathematical details of the 
solution of the eigenvalue equation (4) with adjacency matrix (5); here, we only describe the 
properties of the spectrum which turn out to be relevant for our subsequent analysis. 

The spectrum a is formed by N s states and is divided in three parts: E , a and E + . 
do describes the delocalized states with energies between —It and 2t; in the thermodynamic 
limit L — > oo, the normalized density of these states is given by 
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just as for a particle hopping on a linear chain. Apart from the (pL — 1) continuum states 
belonging to <7o, there are two bound states confined away from the continuum corresponding 
to energies E < —2t and E + > 2t. These two eigenstates are localized and form the so-called 
hidden spectrum [2,3,16]: hidden means here that the two states - in the thermodynamic 
limit - do not contribute to the normalized density of states yielding the closure relation. 

Eq is the ground-state energy which, in the thermodynamic limit and for a star graph 
with p arms, is given by 



Equation (7) reproduces exactly the known result E = —2t for p = 2 and it implies the well 
known fact that, for a linear chain, there are no localized states. The energy of the isolated 
eigenstate in the high-energy region is simply given by E + =\ E \. 

It is worth observing that the spectrum is gapped: in fact, there is a finite gap A = 
l-Eol ~ 2£ between the ground-state energy and the continuum part of the spectrum; the 
value of A depends on the number of arms p. In Fig. 2 we plot the energy gap as a function 
of the number of arms. As expected, one has A = when p — 2. 

The ground-state and the eigenstate corresponding to the eigenvalue E + - due to the 
topology of the array - are localized in the center and exhibit an exponential decay in the di- 
rection of the arms. In the thermodynamic limit, the normalized ground-state wavefunction, 
as a function of the number of the arms p, is given by: 



(the normalization is, of course, over the whole network: i.e., J2 x , y \i>E {x, y)\ 2 = 1)- £ 
provides an estimate of the ground-state localization and it is given by 




(8) 
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(9) 



log (p - 1) " 
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In Fig. 3 we plot the ground-state wavefunction for different values of the number of arms p. 
Figure 3 evidences that adding arms enhances the localization of the wavefunction around 
the center of the star graph. As we shall see in the next section, bosons are allowed to 
spatially condense in this ground-state at low temperatures. 



III. THERMODYNAMICS OF BOSONS HOPPING ON A STAR GRAPH 

The thermodynamical properties of non-interacting bosons hopping on a star-graph hint 
to the possibility of a topology induced spatial BEC in the center of the star graph [3]. To 
elucidate this phenomenon, it is most convenient to introduce the macrocanonical ensemble 
to determine the fugacity z as a function of the temperature of the system. The equation 
determining z is given by 

* - E z -Ji E i.y <«>) 

E £<T 

In Eq. (10) d(E) is the degeneracy of each single-particle eigenstate, E is the energy of the 
ground-state of the Hamiltonian (3) and j3 = 1/JcbT; the sum in Eq. (10) is over the entire 
spectrum a. For free bosons hopping on a star graph, one has 

N T = N Ea + N Et+ J Ee JE z _^ E l_ r (11) 

where N Eo (p, L;T) and N E+ (p, L;T) denote, respectively, the number of bosons which oc- 
cupy the ground-state and the state of energy E + at a certain temperature T. p(E), with 
E G (To, is the density of delocalized states defined in Eq. (6). It is pertinent to de- 
fine also the number of particles per site in each part of the spectrum as = N Eo /N s , 
n ao = J Eeao dE p(E) [ z -i e P&-Eo) - and n E+ = N E+ /N S . In the thermodynamic limit, 
one has 

n Eo (T) = lim 



L^oo N s Z- 1 - 1 

and n E — since 
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n E+ {T) = lim 



l-oo N s Z ~ l e- 2 P E ° - 1 e- 2 ? E ° - 1 








VT. 



Thus, in the thermodynamic limit, E + is not macroscopically occupied at any temperature 
and does not play any role in describing the thermodynamics of the system. 

The last term of the right-hand side of Eq. (11) represents the number of bosons in the 
delocalized (chain-like) states. The presence of the hidden spectrum changes the behavior of 
the integral evaluated in the interval {— 2t, 2t}, since it reduces it to the one describing non- 
interacting bosons on a linear chain with an impurity in one of the sites. As a result, letting 
z — > 1, the integral converges, even at finite temperatures, making possible the topology 
induced spatial BEC in the center of the star graph. 



If one defines T c as the critical temperature at which BEC occurs, for any T < T c , 
the ground-state is macroscopically filled. Since, at the critical temperature and in the 
thermodynamic limit, n Eo (T c ) = 0, from Eqs. (6) and (11) one has that the equation 
allowing to determine T c as a function of the filling / and of the hopping strength t reads as 



Equation (12) can be solved numerically for any value of /. When / ^> 1, one may expand 
the exponential in Eq. (12) to the first order in the inverse of the critical temperature T c 
getting 



A. Critical temperature and condensate fraction 




(12) 




(13) 



Substituting cos6> = E/2t in Eq. (13), one has 




from which 
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k B T c = Ej 



with Ej the Josephson energy defined in Eq. (2). The result (14) holds for any graph for 
which Eq < —It and the density of states of the continuum part of the spectrum is given by 
Eq. (6), i.e., the density of states of a linear chain. Since, for a comb lattice E = — 2y / 2t 
[2], for this graph one gets k B T c = Ej. For a linear chain one has instead E = —It and 
thus T c = 0: of course, no condensation occurs in this case. 

Upon inserting the value of the ground-state energy (7) in Eq. (13), the critical temper- 
ature T c is given by 

kB% * 270i Ej - (15) 

Equation (15) has been checked numerically and it is in excellent agreement with the nu- 
merical solution of Eq. (12): for / 3> 1, the error is of order 1/ /. For interwell barriers V 
of order ~ 15.Br ~ 2n% ■ 50kHz and for fillings / ~ 200 one has Ej ~ 50nK. According to 
Eq. (15), one then expects the formation of an observable condensate in the star center. In 
Fig. 4 we plot the critical temperature T c given by Eq. (15) as a function of the number of 
arms. 

One may use Eq. (15) to determine also the condensate fraction as a function of the scaled 
temperature T/T c . In the thermodynamic limit, the number of particles in the delocalized 
states is given by 

r 2t dE T 

N n (T/T c ) = lim N, j_ 2i p(E) emE _ E(>) _ i „ N T . -. (16) 

In Eq. (16) the exponential has been expanded to the first order in (3: this approximation 
holds for / ^> 1 and it is in very good agreement with the numerical evaluation of the 
integral (16) also in a large neighborhood below T c . The critical temperature at which BEC 
occurs crucially depends on the number of arms of the star (see Eq. (15)) and thus one may 
adjust it by choosing a pertinent number of arms. 

From Eqs. (11) and (16), one gets the number of particles in the localized ground-state 
N Eo : the fraction of condensate, for T < T c , is then given by 
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For / ranging from 10 3 to 10 9 , the results provided by Eq. (17) differ from those obtained by 
the numerical evaluation of N Bo from Eq. (11) by less than 1%. Equation (17) clearly shows 
that the condensate has dimension 1; cigar-shaped one-dimensional atomic Bose condensates 
support, in fact, a condensate fraction given by Eq. (17) [17,18]. 

B. Distribution of bosons along the arms of the star network 

In the following we shall determine the distribution of the bosons over the star graph. 
Due to the topology induced spatial condensation in the center of the star graph, one should 
expect an inhomogeneous distribution of the bosons along the arms of the network. The 
average number of bosons N B at a site (x, y) depends - due to the symmetry of the graph - 
only on the distance x from the center of the star. At any temperature, N B is given by: 

, rlt 1 . 

N B (x;T/T c ) = Km {n Eo (x;T/T c ) |^ (x)| 2 + (pL + 1) dE p{E) -^-^ - { |^)| 2 }. (18) 

In Eq. (18) ipE (x) is the wavefunction corresponding to the ground-state of the single- 
particle spectrum and iPe(%) is the wavefunction associated to a delocalized state with 
energy E. The last term in the right-hand side of Eq. (18) gives then the contribution 
coming from the delocalized states, which, in the thermodynamic limit, is independent from 
the site index x and equals the constant (N T /N S ) ■ (T/T c ). Using Eq. (17), for T < T c , one 
has 

N B Tin „ N T {(l - 1) £^e-*-» + . |} . ( 19) 

The exponential behavior of the ground-state eigenfunction leads, for T < T c , to an increase 
of N B on the sites near the center of the star while, when x ^> 1, the behavior is dominated 
by the last term in the right-hand side of Eq. (19). Thus, away from the center, once the 
filling is fixed, N B depends only on the scaled temperature T/T c and it is given by 
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Topology induced spatial BEC in a system of non-interacting bosons hopping on a star 
graph predicts then a rather sharp decrease of the number of bosons at sites located away 
from the center. The linear dependence exhibited by the solid line in Fig. 5 is consistent 
with the observation that, in this system, the condensate has dimension 1. 



IV. CONCLUDING REMARKS 

We showed how the topology of an optical lattice confining a dilute atomic gas may 
catalyze the existence of new and unexpected coherent phases. For this purpose we analyzed 
the paradigmatic and simple example of bosons hopping on star graph; our analysis allowed 
not only for the computation of the critical temperature T c for which there is condensation 
of the bosons in the center of the star, but also allowed us to compute the distribution of the 
bosons along the arms of the star and to show the simple dependence of the non-condensate 
fraction on the reduced temperature T/T c . We find T c oc tf where t is the tunneling rate 
and / the average number of particles per site. 

In this paper we analyzed the behavior of bosons hopping on a star shaped network in 
the thermodynamic limit. It is comforting to observe that numerical simulations point out 
to the fact that - already for / ~ 100 and for a reasonable number (L ~ 50) of lattice 
sites on each arm of the star graph - finite size effects are negligible and, thus, the results 
derived in this paper are also very useful to pertinently describe the variety of experimentally 
accessible systems for which it is expected to observe the signature of a topology induced 
BEC. Although we focused our attention to bosons hopping on star shaped optical networks, 
the reader may easily convince her(him)-self that our analysis could be also applied to the 
description of topology induced coherent phenomena in star shaped Josephson junction 
networks [2]. In the latter application, there is practically no limitation on the number of 
junctions needed to build the star shaped network and, thus, the results obtained in the 
thermodynamic limit are expected to be very accurate. 
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An experimental realization of a star-shaped optical lattice with four (or six) arms may be 
achieved by first creating a regular square (or cubic) lattice using two pairs of counterprop- 
agating laser beams; for these configurations, the optical potential has the form V(x,y) = 
Vo[sin 2 {kx) + sin 2 (ky)] (or similarly V(x, y, z) = Vo[sin 2 (kx) + sin 2 (ky) + sin 2 (kz)]). A row 
and a column of the square lattice (or three perpendicular chain of the cubic lattice) may 
then be selected by superimposing two (or three) perpendicular gaussian laser beams, ob- 
taining by this procedure a four- (six-) arm star optical lattice. Typical experimental values 
of Vo for which the tight-binding approximation and the Bose-Hubbard Hamiltonian (1) are 
valid are Vo ~ 10 — 30Er (Er = h 2 /2m\ 2 with A the wavevector of the lattice beams, with 
A ~ 800nm). With an average number of particle per site / ~ 200, from Eq. (15) one finds 
(for Vo ~ 15E E and p = 6) T c of the order of 50nK. Below T c , the macroscopic occupation of 
the ground-state could be evidenced by turning off the magnetic+optical trap and observing 
the gas expansion, as in the usual detection of atomic Bose-Einstein condensates. 

The experimental observation of topology induced spatial BEC on a star shaped network 
may be easily achieved also using superconducting Josephson junctions. For a superconduct- 
ing network, it is sufficient to measure the I- V characteristic of a single arm of the Josephson 
junction network (JJN) built on a star graph and of the Josephson critical current along a 
given arm; if, in fact, one feeds an external current I ex t at the extremities of the arm, one 
expects to observe no voltage unless I ext is larger than the smallest of the critical currents 
of the junctions along the arm. Since, below T c , the Josephson critical current of the arm is 
given by the smallest of the critical currents of the junctions positioned along the arm, the 
measurement of the I-V characteristic of an arm of the star graph should provide a measure- 
ment of the critical current of the junction which is farther from the star center. From the 
analysis carried out in this paper, it is rather easy to provide an estimate of the Josephson 
critical current as a function of both the temperature and the distance from the star center 
above and below T c . Furthermore, one may show that the ratio of the Josephson critical 
currents of a Josephson junction - located on a given arm between the sites (x + l,y) and 
(x, y) - above and below the critical temperature T c does not depend on y and it is given by 
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If(x,r)/I^(x) ~ \J N B (x + 1; t)N b (x; t)/ f, where If (I^) is the Josephson critical current 
below (above) T c , r = T/T c is the scaled temperature and Nb(x;t) is given by Eq. (20). 
Far away from the star center (x >> 1), one gets 

Mx) (21) 
Thus, BEC in a star shaped JJN predicts a sharp decrease of the Josephson critical current 
for a junction located away from the center. 

The striking and intriguing similarities between superconducting Josephson junction net- 
works and atomic gas in suitable deep optical lattices have been already pointed out [6]: one 
may think, in fact, to realize a star shaped network also using bosonic Josephson junctions. 
For this purpose, it is needed that, in each well of the periodic potential, there is a con- 
densate grain appearing at a Bose-Einstein condensation temperature Tbec an d that, when 
Tbec is larger than all other energy scales, the atoms in the i-th. well of the optical lattice 
may be described by a macroscopic wavefunction fy. Thus, it becomes apparent that an 
optical network may be regarded as a network of bosonic Josephson junctions. Our analysis 
shows that - at a temperature T c < Tbec - the topology of the star-shaped network induces 
a further finite temperature transition to a state in which the bosons spatially condense in 
the center of the star. 
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APPENDIX A: ENERGY SPECTRUM ON A STAR-SHAPED NETWORK 

In this Appendix we solve the eigenvalue equation (4) with the adjacency matrix given 
by Eq. (5), describing bosons hopping on a star graph with (pL + 1) sites (i.e., a star with 
p arms having L sites each). 
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Since non-interacting bosons on a linear chain are described by plane waves with wave 
vector k one has that, on each arm, an eigenstate of Eq. (4), corresponding to energy 
E = —2tcos(k), may be written as: 

i/>(x, y) = A y e tkx + B y e~ ikx . (Al) 

In Eq. (Al) y — 1, • • • ,p is an index labeling the arm, while x — 0, • • • , L labels the sites on 
the arm. The wavefunctions described in Eq. (Al) are, of course, delocalized. 

Requiring that - on each arm - ip(x, y) is a solution of the eigenvalue equation at x = L 
amounts to require that A y , B y and k should satisfy the p equations 

-t (A y e lk{L - l) + B y e^ L -^) = -2t cos(k) (A y e ikL + B y e~ ikL ) y = l,---,p. (A2) 

Furthermore, the eigenstates defined in Eq. (Al) should satisfy (p— 1) matching conditions 
in the center of the star where the wavefunctions defined on each arm are linked; thus, one 
has 

A y + B y = A y+1 + B y+1 y = l,---,p-l. (A3) 
Using Eqs.(A3), the condition in the center gives one more equation 

-t^(A y e ik + B y e~ lk ) = -2t cos(k) (A y , + B y ,) y' = l,---,p. (A4) 

y=l 

Equations (A2) and (A4) may be grouped in a homogeneous linear system of 2p equations 
which allows to fix the 2p parameters A y and B y . Upon denoting with M the (2p x 2p) 
matrix whose elements are the coefficients of the linear system given by Eqs. (A2) and (A4), 
requiring that 

det M = Q{k, L) ■ (1 - e ™(L+i)y-i .{(p-2) cot{k) - pcot [k{L + !)}} = 0, (A5) 

guarantees the uniqueness of the solution. In Eq. (A5) \Q(k, L) \ = 1 for any value of k. 

One immediately sees that the values of k for which k = nn/(L + 1) (with n = 1,2, L) 
provide a set of L {p — l)-fold degenerate eigenstates of Eq. (A5). In addition, the solutions 
of the transcendental equation 
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(p-2)cot(fc) -pcot[k(L + l)] =0 (A6) 

provide the values of k associated to non-degenerate eigenstates. Equation (A6) can be 
solved numerically and yields a set of (L — 1) non-degenerate eigenvalues corresponding to 
values of k which - in the thermodynamic limit - are equally spaced and separated by a 
distance n/(L + 1). As a result, the set of delocalized states is formed by (pL — 1) states 
corresponding to energies ranging between —It and +2t. One can easily convince oneself [3] 
that, in the thermodynamic limit, the normalized density of states is given by 

" (E) = WW=W' (A7) 

as in the case of non-interacting bosons hopping on a linear chain (see e.g. [13]). 

Since the total number of states should equal N s = (pL + 1), there are also two localized 
states in the spectrum: to find them, it is convenient to look for solutions of the eigenvalue 
equation (4) of the form 

ifo(x) = Ae- vx + Be vx 

if}+(x) = A(-l) x e-^ + B(-l) x e vx (A8) 

corresponding, respectively, to the eigenvalues E 0>+ = =F2tcosh?7. In Eqs. (A8) A and B 
are normalization constants and r\ = l/£ is a parameter accounting for the localization of 
the states. One may determine the parameters A, B and r\ by using only the normalization 
condition for ifj and by rewriting Eq. (A6) for k = irj. Namely, one should solve: 

(p - 2) coth(?7) - p coth[??(L + 1)] = 0, (A9) 

together with the condition that J2 x , y l^ot^? y)\ 2 — 1- For L — > oo, one can always set B = 0. 
Equation (A9) becomes (p — 2) coth(^) — p = which is solved by r\ = log(p — l)/2 yielding 
E = —tpj y/p — 1 and E + = —E . Solving the eigenvalue equation (4) for E = E 0<+ , one 
obtains the wavefunctions of both the localized states. The normalized eigenfunction for the 
ground-state is then given by: 
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^ x ) = {i^ e ~ x/ *- ( A1 °) 

Since, for L — > oo, the normalized density of the continuum states is given by p(E) = 
l/(ny/4:t 2 — E 2 ), the two localized states do not contribute to the closure relation in the 
thermodynamic limit; thus, they belong to the hidden spectrum. 
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FIGURES 




FIG. 1. A star network with six arms. 
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FIG. 2. Energy gap A between the ground-state energy and the continuum part of the spectrum 
(in units of the tunneling rate t) as a function of the number of arms p. For p = 2 (i.e., a linear 
chain) A = 0. The dotted line is just a guide to the eye to connect the points. 
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FIG. 3. The normalized single-particle ground-state wavefunction for bosons hopping on a star 
graph as a function of the distance x from the center. The number of arms p is respectively 3 (solid 
line), 6 (dotted line), and 10 (dashed line). 




°2 ' 4 ' 6 ' 8 ' 10 

P 

FIG. 4. Critical temperature T c (in units of ks/Ej) as a function of the number of arms p. For 
the linear chain (p = 2) T c = 0. 
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FIG. 5. Distribution of the number of bosons Nb as a function of T/T c computed for x » 1. 
Nb(x) is in units of the filling / and is therefore equal to 1 for T > T c . 
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